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Abstract 

An important class of physical systems that are of interest in practice are input- 
output open quantum systems that can be described by quantum stochastic dif¬ 
ferential equations and defined on an infinite-dimensional underlying Hilbert space. 

Most commonly, these systems involve coupling to a quantum harmonic oscillator 
as a system component. This paper is concerned with error bounds in the hnite- 
dimensional approximations of input-output open quantum systems defined on an 
infinite-dimensional Hilbert space. We develop a framework for developing error 
bounds between the time evolution of the state of a class of inhnite-dimensional 
quantum systems and its approximation on a finite-dimensional subspace of the orig¬ 
inal, when both are initialized in the latter subspace. This framework is then applied 
to two approaches for obtaining finite-dimensional approximations: subspace trunca¬ 
tion and adiabatic elimination. Applications of the bounds to some physical examples 
drawn from the literature are provided to illustrate our results. 

Keywords-. Quantum stochastic differential equations, input-output open quantum systems, 
finite-dimensional approximations, error bounds, approximation errors 


1 Introduction 

Quantum stochastic differential equations (QSDEs) developed independently by Hudson 
and Parthasarathy [1] and Gardiner and Collett [2] (the latter in a less general form than 
the former) have been widely used to describe the input-output models of physical open 
markov quantum systems HHS]. Such models describe the evolution of Markovian quantum 
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systems interacting with a propagating quantum field, such as a quantum optical field, 
and are frequently encountered in quantum optics, optomechanics, and related fields. An 
example in quantum optics would be a cavity QED (quantum electrodynamics) system 
where a single atom is trapped inside an optical cavity that interacts with an external 
coherent laser beam impinging on the optical cavity. These input-output models have 
subsequently played an important role in the modern development of quantum filtering 
and quantum feedback control theory m. Many types of quantum feedback controllers 
have been proposed in the literature on the basis of QSDEs, using both measurement- 
based quantum feedback control, e.g., [3|[T1[6|[7], and coherent feedback control, e.g., [HHin]. 
Besides, QSDEs have also been applied in various developments in quantum information 
processing, such as in quantum computation technology; e.g., see ra¬ 
in various physical systems of interest, one often deals with input-output systems that 
include coupling to a quantum harmonic oscillator. For instance, typical superconducting 
circuits that are of interest for quantum information processing consist of artificial two-level 
atoms coupled to a transmission line resonator. The former is typically described using a 
finite-dimensional Hilbert space and the latter is a quantum harmonic oscillator with an 
infinite-dimensional underlying Hilbert space (i.e., L^(R), the space of square-integrable 
complex-valued functions on the real line). Another example is a proposed photonic realiza¬ 
tion of classical logic based on Kerr nonlinear optical cavities in |T2], which is built around 
a quantum harmonic oscillator with a Kerr nonlinear medium inside it. If a mathemati¬ 
cal model for such quantum devices is sufficiently simple, it is often possible to simulate 
the dynamics of the system on a digital computer to assess the predicted performance of 
the actual device, as carried out in [12]. The simulation carried out is typically that of a 
stochastic master equation that simulates the stochastic dynamics of a quantum system 
when one of its output is observed via laboratory procedures such as homodyne detection 
or photon counting, see [SlIEllIS]. However, since it is not possible to faithfully simulate a 
quantum system with an infinite-dimensional Hilbert space, often in simulations this space 
is truncated to some finite-dimensional subspace. Two approaches that are often employed 
to approximate a quantum system, defined on an infinite-dimensional space, are subspace 
truncation approximation and adiabatic elimination (also known as singular perturbation). 
Subspace truncation approximation is applied to eliminate higher dimensions of the origi¬ 
nal infinite-dimensional Hilbert space. An operator X on the infinite-dimensional space is 
approximated by a truncated operator of the form PXP, where P denotes an orthogonal 
projection projector onto the approximate finite-dimensional subspace. For instance, with 
quantum harmonic oscillators, a commonly used finite-dimensional space is the span of a 
finite number of Fock states |0), |1),..., |n). On the other hand, adiabatic elimination is 
often used to simplify quantum systems comprising components that evolve at multiple 
well-separated time-scales. In this approach, the faster variables are eliminated from the 
mathematical model description of the systems. 

Despite the ubiquity of approximating infinite-dimensional Hilbert spaces of quantum 
systems by finite-dimensional subspaces for simulations of input-output quantum systems, 
to the best of the authors’ knowledge, there does not appear to be any work that has tried to 
obtain some explicit bounds on the approximation error of the joint state of the system and 
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the quantum field it is coupled to. This work develops a framework for developing bounds 
on the error between the quantum state of a quantum system described by the QSDE and 
the quantum state of a finite-dimensional approximation described by another QSDE, when 
both systems are initialized in a state in the finite-dimensional subspace. Central to the 
framework is a contractive semigroup associated with the unitary QSDEs of input-output 
Markov quantum systems. Error bounds are developed for both adiabatic elimination 
and subspace truncation approximations. For illustration, our results are applied to some 
physical examples drawn from the literature. Prelimary results of this work were announced 
in the conference paper na. The results presented in the work go significantly beyond na. 
In particular, na only treates the subspace truncation approximation with some elements 
of the proofs omitted, error bounds for adiabatic elimination had not been developed, and 
computability of the error bounds were not considered. 

The rest of this paper is structured as follows. In Section [2l we present the class of 
open quantum systems and the associated QSDEs describing Markovian open quantum 
systems. Explicit error bounds for the subspace truncation approximation of a Markovian 
open quantum system are established in Section [3] and some examples are provided. We 
then establish error bounds for adiabatic elimination approximation in Section 0] and some 
examples are also provided. Finally, concluding remarks close the paper in Section [5l 


2 Preliminaries 

2.1 Notation 

We use i = \/^ and let (•)* denote the adjoint of a linear operator on a Hilbert space as well 
as the conjugate of a complex number, and denote matrix transposition. We denote 
by 6ij the Kronecker delta function. We define = ^{A + A*) and —H*). 

For a linear operator A, we write ker(H) to denote the kernel of A and ran(H) the range 
of A. We often write j-) to denote an element of a Hilbert space and denote by the 

algebraic tensor product of Hilbert spaces Ti and J^. For a subspace "Hq of a Hilbert space 
"H, we write P-Ho to denote the orthogonal projection operator onto T-Lq. For a Hilbert space 
"H = 'Ho © Hi, we will write H 0 Ho to denote Hi. For a linear operator X on H, X\^^ 
denotes the restriction of X to Hq. We use H(H) to denote the algebra of all bounded 
linear operators on l-i. We write [H, B] = AB — BA. The notation || ■ || will be used 
to denote Hilbert space norms and operator norms, (•, •) denotes an inner product on a 
Hilbert space, linear in the right slot and antilinear in the left, and j-) (-1 denotes an outer 
product. Here, l[o,t] : [0,f] —)■ {0,1} denotes the indicator function. Finally, Z_|_ denotes 
the set of all positive integers. 

2.2 Open quantum systems 

Consider a separable Hilbert space Hq and the symmetric boson Fock space (of multi¬ 
plicity m) X defined over the space L^([0, T]; C™) = C™ © L^([0,T]) with 0 < T < cxo; 
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see [ini Ch. 4-5] for more details. We will use e(/) G -F, with / G L^([0, Tj; C™'), to 
denote exponential vectors in F”. Let & C L^([0, Tj; C™) fl L“^([0, Tj; C™') be an admissi¬ 
ble subspace in the sense of Hudson-Parthasarathy [T] which contains at least all simple 
functions, where L“(,([0, Tj; C™) is the space of locally bounded vector-valued functions. 
Here, we will consider a dense domain Pq C 'Hq and a dense domain of exponential vectors 
£ = span{e(/) | / G ©} C F”. 

Consider an open Markov quantum system that can be described by a set of linear 
operators dehned on the Hilbert space (i) a self-adjoint Hamiltonian operator H, (ii) a 
vector of coupling operators L with the j-th element, Lj : Hq —)■ Hq for all j = 1, 2,..., m, 
and (hi) a unitary scattering matrix S with the ij-th element, Sij : Hq —?• Hq for all i,j = 
1,2,... ,m. Moreover, the operators Sij, Lj, H and their adjoints are assumed to have Vq 
as a common invariant dense domain. Under this description, we note that m corresponds 
to the number of external bosonic input helds driving the system. Each bosonic input 
held can be described by annihilation and creation held operators, 6) and bf, respectively, 
which satisfy the commutation relations [6j, bi ] = 6ij6{t — s) for all i,j = 1,2,... ,m and 
all f, s > 0. We can then dehne the annihilation process M), the creation process Al , and 
the gauge process as 

bids, Ai* = [ bi*ds, 

Jo 

Note that these processes are adapted quantum stochastic processes. In the vacuum repre¬ 
sentation, the products of their forward diherentials dAl = Al_^_^^ — Al, dA]* = Al_^_^f* — A]*, 
and dAl^ = — A]^ satisfy the quantum ltd table 


A« = r K'Hds. 
Jo 



X 

dA'l 

dA'p* 

dAf 

dt 

dA\ 

0 


dikdA\ 

0 

dAi* 

0 

0 

0 

0 

dh\’ 

0 

djkdAf 

SjkdAf 

0 

dt 

0 

0 

0 

0 


Here, b] = can be interpreted as a vacuum quantum white noise, while A“ can be 
interpreted as the quantum realization of a Poisson process with zero intensity pQ. 

Following ra. the time evolution of a Markov open quantum system is given by an 
adapted process Ut satisfying the left Hudson-Parthasarathy QSDE [1]: 


dU, = f/d Xi - d<i)dtdf) + E 


.*J = 1 


2=1 






2=1 


dt 


( 1 ) 


with Uq = I. The quantum stochastic integrals are dehned relative to the domain Vq^S. 
With the left QSDE, the evolution of a state vector 'ip G Hq ® H is given by 
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In this paper, we are interested in the problem of approximating the system with 
operator parameters {S, L, H) by an open qnantum system with linear operator parameters 
dehned on a closed subspace C Ho, where is unitary and is 
self-adjoint. Consider a dense domain C H^^\ Again, the operators (S^^\ 

and their adjoints are assumed to have as a common invariant dense domain. Similar 

to o. the time evolution of the approximating system is given by an adapted process 
satisfying the left Hudson-Parthasarathy QSDE |1]: 


{ m m 

^ - %)dA«) + (if’v.4; 

i,i=l i=l 


*J=1 
(k)* 


- E [srL.dAr) + 

*j=i 




i=l 


dt > , 


( 2 ) 


ik') 

with Uq = I. Here, the quantum stochastic integrals in the above equation are dehned 
relative to the domain Similarly, the evolution of a state vector xb G ® is 

given by Ul yj. 


2.3 Associated semigroups 

Let 9t : L^([f, T]; C™') L^([0, T]; C™') be the canonical shift 6tf{s) = f{t + s). We also 

let Qt '■ ^ ^ denote the second quantization of dt, where T\^t denotes the Fock space 

over oo)] C”^). Note that an adapted process Ut on Hq ^ H is called a contraction 

(or unitary) cocycle if Ut is a contraction (or unitary) for all f > 0, f i—)■ Ct is strongly 
continuous, and f/^+t = f/s(/ ® Q*UtQs)- 

Let us now impose an important condition on the open quantum systems under con¬ 
sideration, adopted from [TB] . 

Condition 1 (Contraction cocycle solutions). For all t > 0 and all k G 


(a) the QSDE ([T]) possesses a unique solution Ut which extends to a unitary cocycle on 
Ho ® H, 

(b) the QSDE ([2]) possesses a unique solution Ut^'^ which extends to a contraction cocycle 
on H^^'> (8) H. 

Let us dehne an operator : Ho -> Ho via the identity 

{u, = e"5(ll"lh+ll/3|lbi (g) e{alyo,t]), UtV (g) e(/?l[o,t])) 


for all G Ho and aW a, (3 G From [161 Lemma 1], under Condition [11(a), the 
operator G B(Ho) is a strongly continuous contraction semigroup on Ho and its 
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generator satisfies Dom(£("^)) D Vq such that 


5; (a-S*.ft) - 5; (a’SJij) + 5; (L*ft) 

Jj = l *J = 1 j = l 

^ + iI/9|K 






a 


i=l 


U 


( 3 ) 


for all u G Pq- Here, we note that Doni(£’'"^)) is dense in "Ho- We likewise dehne an 
operator —). by replacing Ut with 

Condition 2 (Core for generators). For allt>0 and all k G 

(a) Vo is a core for 

[b) V^^') is a core for 


Condition l^g) ensures that the dehnition (|3]) completely determines £(“^1 for all a, /9 G 


C™, likewise 1 ^6) [ completely determines In the sequel, we will make use of the above 

semigroups associated with open quantum systems in establishing our model approximation 
error bound. Several sufficient conditions are known to guarantee that a QSDE possesses 
a unique solution that extends to a unitary cocycle when the Hilbert space is inhnite- 
dimensional and the operator coefficients of the QSDE are unbounded, see, e.g., na, and a 
related discussion in [Ml Remark 4], Throughout the paper, we will assume that Conditions 
[T] and [2] are fulhlled. 


3 Error bounds for subspace truncation approxima¬ 
tions 


In this section, we consider the problem where the inhnite-dimensional space T-Lq is trun¬ 
cated to a hnite-dimensional subspace and the original operators X on FLo is ap¬ 

proximated by truncated operators of the form Py^(k)XP^(k). Here, the dimension of 
increases with k G Z_|_ and V^^l = Moreover, Condition [1 1(6) (in fact, is uni¬ 

tary [1]) and E 1(6) hold immediately. 


3.1 Assumptions and preliminary results 

Assumption 1. For any fc G and any a,/3 G C™, FL^^'> C Dom(£(“^)). 

Let = ran (_ £(*;;«,/3)^ |-^(fe))- Supposing that Assumption [1] holds, we also 
assume the following. 

Assumption 2. For each fc G and each a, fd e C^, there exists 7^“^^ > 

0, and a non-trivial subspaee {0} C C such that 
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(a) 

(b) (£M) - = 0, 

z.e., •H(^) © /C(^) C ker ((£(“« - 

(c) For any u G 

^icww)} = -9{k, a, /3) ||P^(fc)M||^ + h{k, a, /3, m) 


/or some 5 f(fc, a, /3) > and some \h{k, a, /9, m)| < II-Pac('')'*^II II'^^II ■ 

(d) For any u G and any t > 0, 

S{{T©>£l”'’>P^,„t.,T7’«P^,.,«>} = -?j(k,aJ) 




+ h{t, k, a, 13, u) 


ye 




for some g{k, a, (3) > and some h{t, k, a, (3, u) 

Assumption 3. There exists r,s E Zi+ such that, for all a, (3 E C™, we have that 


\u\\. 


lim I — 

fc->oo 


/Jk;ay)\ 

qa 


(afS) 

Tk 


(a/3) 

% 


= 0 


Moreover, for any i = 0,1,2,... ,r, and j = 0,1,2,..., s, we have that 


Tc 


{k-,alS) / (fc;a/3)\ / (fc;ay) \ F 


lim , 

k^oo 


Qe 


{a/3) 

Tk 


{a/S) 

Tk 


= 0 


{k-,ay) {k-,ay)\ ( 1-2 
(fc;a/3) / <?e qa 


lim e '1* 


/c—>-oo 


( 7 , 


(«/3)\2 


= 0 , 


for any a, (3 E PF, any t > 0, and for i = 0,1,2,..., min{r, s} — 1. 


Let us present some useful lemmas. 


Lemma 1. Suppose that Assumption \ ^c) \ holds. Then for any k,r E Tj+, any a, (3 E C™, 
any u G p(^(, and any t > 0, zt holds that 


PfcwTi 


< 


Jpo.3)' 

qa _ 

(a/3) 

% 


1 - 2 - 


r—1 

j=0 


I fa 


{k;al3)' 


(1-2-0' 


(a/3) 

Tk 


\u\ 


(4) 


where Cq = 1 and Cj 
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Proof. First note, from the definition of a strongly continnous semigronp, that = / 

and [18]. From Assnmption ^|(c)[ we have that 

d 
dt 


p rp(k-,afi)^ 


^ ^ / D rp(k\af}) 73, rp(k-,af3) 

= ^ {PicwTi ^’u 


.^{k-,alS) 


)c{k)lt u 

= 23fJ P,^(k)Tf’-'u 

=-2g{k, a, (3) ||p^wTf V 2h{k, «, /3, Pf 
Solving the above ODE gives ns that 

= \\PKWuf + 2 [ h{k, a, /3, pf 

Jo 

< llPjcwPf||n||dr (5) 

Jo 

{k\afi) 




< 




(o/3) 

7fc 


M 


Here, the second step follows from Assnmption [ ^c)[ The last step follows becanse ||P;s;(fc)M|| < 
and is a contraction semigronp. Noticing that y/laP + |6P < |a| + \b\ for any 


a, 6 G M, we have that 


Pic(k)Tl 

< 

( m f \ 2 

Ma 

^ ^(a/3) I 



[ \7a: / J 


\u\ 


Now, snbstituting (j6|) into the right-handed side of ([5|), we have that 
2 


( 6 ) 




<e-^-^“\uf + 2q, 


(k]OL^) 


e '^^WuW'^dr 


{k]OL^) \ 2 


+ A-5- Wufdr 


(afi) 

Tk 


{k]OL^) ' 


(l - ( ^ ) + (l - 


-2-vF«A / 


< 






7fe 


(q/3) 



(fc;o/3)\ ? 


(«/ 3 ) 

7fc 




Taking the sqnare root on both sides of the eqnation, we get 


P/cwT^ ’u 


< 


(k](y.^) \ 2 


(ckt^) 

7fc 


\u\ 




{k\aj3) \ 2 


7r> 


{k\a0) \ 4 

ya_ 

{af}) 

Tk 


\u\ 














































From repeat application of the above steps, we establish the lemma statement. □ 

Lemma 2. Suppose that Assumption holds. Then for any k,r E , any a,f3 E 
any u E and any t > 0, it holds that 




< 


% 




ye 


1-2-^ 


j=0 


% 


l«ll (7) 


where Cq = 1 and Cj = \JCj-i2i (2i — 1) ^ for j > 1. 


Proof. Similar to Lemma [H using Assumption E (c?), we have that 


2 

= —2g{k,a,ft) Pj^(k)U + 2h{t,k,a, ft,u). 

Solving the ODE, we have that 




„(a/3) 


(°^P) ! 






< 




{k;al3) 


{a/3) 

Tk 


\u\ 


The lemma statement is then established by following similar arguments to Lemma [TJ □ 


3.2 Error bounds for finite-dimensional approximations 

We begin by dehning 


Jpo-h) 

nc 


Jpo-h)' 

ye_ 

(a/3) 

Tk 


(1-2-0 / (A:;a/3)\ (1-2-0 
qa 


Tk 


(a/3) 


r—1 


2‘ci 

— TT) 

i=0 Ik 




l-e-2 


(fc;o/3)' 


(o/3) 

% 


2T 


S—1 

^ (a^) 

i=l Ik 




l_e-2 'Ik 


(fc;o/3)' 


(i-2-b 


( 1 - 2 -* 


(a/3) 

Tk 



( 1 - 2-0 


(1-2-0 


r—1 s—1 

EE 


2(*+0 


CiCj - 


(“/3) 


.= 0 ,-. (2‘-2J)7/' 

min{r,4-l / (fe;a/3) (fe;a/3) \ (1-2-*) 

..2„-2-Ofc“^h ( T Qa 


Ak;al3) 


\(l-2-0 


,-2“0i“^d _ -2-p^°‘^h\ I Qe \ I qa 


(fc;Q;/3) 


.(1-2-0 


(a/3) 

7fc 


7fc 


(a/3) 


‘ E 


c,e 


i=0 


(7. 


(“/3)^2 


( 8 ) 
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where cq = 1, and cj = Cj-i2^ (2^ — 1) ^ for j > 1. We now establish an error bonnd 

between the two semigronps associated with the open qnantnm systems. 

Lemma 3. Suppose As sumptions \J\ and\E hold. Then for any k,r,s G any a, /5 G C™, 
any u G and any t > 0, it holds that 




(o/3) 


-z 


{k-,afi) 


U 




(9) 


Proof. First note, from the definition of a strongly continnons semigronp, that [IB] 


1. = I for all k e Z+, 

2. for all u G Dom(£(“^)), 

3. _ £{k;a0)y^ fQj. gi^ii g since HZ is finite-dimensional). 

From the above properties and Assnmption [1], we can write for all u G and all t > 0 
that 

— _ rp(k-,al3f^ ^ _ ^£^{ay) jn(a/3) ^(^k-,ap)jn{k-,al3)'^ ^ 

= j y ^ ('^M) _ £(fc;a/3)^ (10) 


with j M = 0. Note that by Assnmption [2], is a 

bonnded operator. Since is a generator of a semigronp on a Hilbert space, 


(^Z{<^p) _ ^(k-,ay)'^ ^ ^ Dom(£(“^)), 

(dne to Assnmption [T|), and 

1^^^^ j.ik;ay)^ ^ C'([O, t]; T^o) 

(the class of continnonsly differentiable fnnctions from [0,f] to PLo), a nniqne solntion of 
ffTOj) exists and is given by [18], Thm 3.1.3] 


^(ap) _ jn(k-,a/3)'^ y _ J rpiap) ^^(a/3) _ £^(k-,ay)^ 

for all f > 0 and all u G PL^. From Assnmption 2 j(6) and the definition of we then 

have for all u G PL^ and alH > 0 that 


n(a/3) rp{k-,a.y) 


U 


< 


^£( o / 3 ) _ £{k-,a0)^ J'P\ay)^ 


dr 




u 


dr. 


( 11 ) 
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Now using the bounds (jl]) and ([7]) (established in Lemmas [T] and [21 respectively), and 


applying Assumption [ ^g) we have that 


< q'c 


{k-,af3) 


He 

{aj3) 

Tk 


( 1 - 2 -’-) 


Qa 


(fc;Q;/3) 


(l-2-'>) 


(ajS) 

% 


(ja 


(1-2-’-) 



r—1 / (fe;a/ 3 )\d 2 )■ 


i=0 

■ s-1 


2=1 


C,;e 




(a/3) 

7 k 

(k;a0)' 


( 1 - 2 - 


(a/3) 

7fc 


r—1 s—1 




qa 


_ i=0 j=0 


{a/3) 

7 k 


{a/3) 

7k 


n • 


The result ([2]) then follows from substitution of the above identity into fllip and integration. 
This establishes the lemma statement. □ 


Let & C L^([0, T]; C™) denote the dense set of all simple functions in L^([0, T]; C™). 
That is, for any f G [0,T] and / G ©', there exists 0 < £ < cxo and a sequence 0 = to < 
ti < ■ ■ ■ < ti < ti+i = t such that / = X]i=osome constants a{i) G C"*, 
i = 0,1,..., £. Let = {u (8) e(/) | u G f G ©'}. We can now proceed to derive 
error bounds for approximations by subspace truncation. 

Lemma 4. Suppose Assumptions Q and H hold. For any "01 = ni 0 e(/i),'02 = ^2 0 
e(/ 2 ) G let to = 0 < ti <...< ti < ti+i = t with 0 < t < T be a sequence 

such that fj = X]i=o Q^j(t)l[u,ti+i) for all i = 0,1,i, with aj{i) G C”^ for j = 1,2 and 
i = 0,1,..., Then for any k,r,s E Z+, 





< '^4,s{ifi+i-ti),ai{i),a2{i)) ||0i||||02 

i=0 


( 12 ) 


Proof. First recall that our admissible subspace & contains &. Hence, C 
and the quantum stochastic integrals are well dehned for all 0 G Also, recall that 

I10IP = {u® e{f),u® e{f)) = ||'n|P||e(/)|0 for any 0 G "Ho <8)-P [I5]. Using the cocycle 
properties (Condition [1]) as well as the dehnitions of and we have the identity 

m 


(t/;0i,02) 


and likewise when Ut and 


= l|e(/,)ll l|e(/ 2 )|| 

are respectively replaced by and 


(13) 
from which 
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we immediately obtain 


Ut - m 




< llV'ill l|e (/2 
Now, note that for any u G that 


ryn(ai (0)02(0)) AT-i(oi(^)a2(r)) 

-‘-tl-to '''-‘-t-ti 


ry-,(fc;oi (0)02(0)) rp{k-,ai{£)a 2 {e))\ 

-‘■tl-to '"-‘t-ti ) 


U2 


(14) 


ry-,(oi (0)02 ( 0 ))ryn(oi (1)02(1)) rp(ai{l)a2{(.)) 

-‘tl-to -‘t2-tl '''-‘t-ti 


ry-,(fc;oi (0)02 ( 0 ))ry-,(A;;oi (1)02(1)) ^(fc;oi (£)a2 (£)) 
-‘tl-to -‘t2-tl "'-‘t-ti 


U 




I 

-r -tf ,— 


(01(0)02(0)) _ ry-i(fc;oi( 0 )o 2 ( 0 ))\ rT-,(fc;0l(l)02(l))rT-|(fc;0l(2)02(2)) ry-i(fc;oi(^)o 2 (f)) 

h-to -‘tl-to I -‘t 2 -tl ^t3-t2 '"-‘t-ti 

(fc;oi(^)o2(^)) 

t2—t\ I t2—t\ t—t£ 


(01(0)02(0)) z'2-1(01 (1)02(1)) _ r7-i(fc;oi(1)02(1))^ r7-i(fc;oi(1)02(1)) 


tl-to 


t 2 -tl 


Ti 


■■■Ti 


I ryn(oi (0)02(0)) rp{oil{e.-l)a 2 {£-t)) (rp(oil{e.)a 2 {£)) _ rp(k\ai(t.)a 2 (t.))\ 

T ^-‘tl-to "'-‘ti-ti_i y^t-ti -‘t-te J 


u. 


From the bound ([9]) (established in Lemma [3]), and that fact that the semigroups are 
contractions, we have for any u G that 


^(01 (0)02 ( 0 ))ry-,(oi (1)02(1)) ry-,( 0 l(£) 02 (^)) 

-‘tl-to -‘t2-tl '"-‘t-ti 


< ^ ((ti +1 - ti),ai{i),a 2 {i)) ||m| 

j=0 


^(fc;oi (0)02 ( 0 ))ryn(fc;oi (1)02(1)) rp{k-,ai{tja2{i)) 

-‘tl-to -‘t 2 -tl "'-‘t-ti 


U 


(15) 


The bound flT^ then follows by substituting flTHl) into flTTl) . This establishes the theorem 
statement. □ 

Corollary 1. Suppose AssumptionsUi IM andlM hold. For any t G [0,T] with 0 < T < 00 , 
any tjji = Ui ® e(/i), -02 = U 2 ® e(/ 2 ) G 0 , we have that 


Ut - ^ 1 ,^ 2 ) < 2(||Mi||||e(/i) - e(/()||||i/> 2 || + ||M 2 ||||e(/ 2 ) - eif^) 

e 


i =0 


for any '0' = uj 0 e(/j) G with /j = fpiMuM+i) some i G Z+, some sequence 
to = 0 < ti <...< ti < ti+i = t, and some constants /L g C™ for j = 1,2 and 
i = 0, l,...,^. Moreover, for any fixed positive integer p G Z+, and any tpi = 0 
e{fi),fj2 = U2® e(/2) G 0 J^, 


lim 

fc—)-oo 


Ut-Ut 


(k) 




= 0 . 


(17) 


Proof. Recall that Ut is unitary and is a contraction (Condition [T]). From the triangle 
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inequality and Cauchy-Schwarz’s inequality, we note that 


U,-U, 


(k) 




< 


Ut - Ml 0 (e(/i) - e{f [)), U 2 0 e(/ 2 )^ 


+ 

(^(Ut - Ml 0 e(/(),M2 0 (e(/2) 

-e(/2))) 

+ 

Ml 0 e(/(),M2 0 e(/ 2 )^ 



< 2||mi|| ||e(/i) - e{f[)\\ \\i) 2 \\ + 2 ||m 2 || ||e(/ 2 ) - e(/2)|| HV'il 

Ut - ui’^A Ml 0 e(/ 0 , M2 0 e(/2) 


The result flTBD then follows from the bound flT^ (established in Lemma 0]). 

To show cm. recall that & is dense in L^([0, T]; C™). Therefore, for each e > 0 and 
each fj G T^([0, T]; C™'), there exists /j G 6' such that ||/j —/j|| < e. Moreover, there exists 

0 < £ < cxo and 0 = to < ^i < • • • < < ti+i = t, such that f' = Yl\=o^'j{^)'^[UM+i) with 

a'(0),a'(l),... ,aj{i) G C™. Suppose that mi,M 2 ^ 0 (otherwise the corollary statement 
becomes trivial), then for any e > 0 we may choose /(,/2 G & (choosing f[ hrst followed 
by f^) such that 

l|e(/i)-e(/()" 


< 


6||mi|| ||e(/2)|| 


|e(/ 2 )-e(/')||< 


6 ||m 2 || I|e(/()|| 

Finally, from Assumption [3] and the bound fll2p (established in Lemma 0]), we can find a 
sufficiently large k G larger than p, such that 




< 


i =0 


3||V'illll^/'2l 


From (061) and the above choices, we then have that 


Ut-U 


(k) 


^ 1,^2 


< e. 


This establishes the corollary statement. □ 

Theorem 1. Suppose that Assumptions [7p|. and 0 hold. Let 0 < T < oo. For any 
t G [0,T], consider any G any ® ^ 

^ 7 ^ 0. Also, consider any f G &. Let i be a positive integer and to = 0 < ti < ... < 

ti < f£+i = t be a sequence such that f = ELo dyt = ELo for 

some constants ff gj^t,i G C™ for for j = 1, 2,..., and i = 0,1,..., i. Let u G with 

||m|| = 1, and I/) = e(/)/||e(/)|| G IF (i.e., |/) is a coherent state with amplitude f), and 
ip = M 0 I/). Then, 

2 


Ut-U 


(fc) 


ft 


L't I 

+ ^r,s {ifi+i - ti), fl dE*) HE,til- 

j=l i=l 


( 18 ) 
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If is unitary for each t then the following bound holds, 


Ut-U, 


(k) 




<4(iii/>-i/'>ii+ c'f’V-v-; 

L't i 

+ 2X] ((^*+1 - ^0, fla'ju) 

j=i i=i 


Moreover, for any fixed p E and any ^jJ = u ® \ f) E 0 J^, 


(19) 


lim 

k^oo 


Ut-U, 


(k) 




= 0 . 


( 20 ) 


for any t E [0, T] with 0 < T < oo. 


Remark 1. Note that a stronger result of strong convergence uniformly over compact time 


intervals, linifc^oo supo<f< 2 ’ 


r(^) 




= 0, has been established in m Proposition 


Ut - m 

20] employing a Trotter-Kato theorem. However, no error bound as in ffTSjl for a finite 
value of k has previously been established. 

Proof. First note that since Ut is unitary and a contraction, 


Ut-Ut 


(k) 


<2m,-uPYi’,um. 


Also, we have that 


+ 


Ui-ulNi’Wli’) < (-0’)’'“»(!/)-I/'))v;v') 

Ut-u'py u®\f'),i.u;Y--iii'S) 

<2|l«l||l/)-|/')IIIIV'll + 2||«|||||/')|| \\U:Y-i>',\ 

p 

+ E (('^•-0’) «®i/'),0) ■ 

i=i 

The result flTK]) then follows from the bound flT^ . and substituting ||u|| = 1, 

')|| = 1. If Ut^'^ is unitary then we have the bound. 


= 1, and 


Ut-Ut 


(k) 




<2K((7,-i7<'-'')-^.,r/f>‘V.)|, 
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and following analogous calculations to the above yields, 


Ut-U, 


(k) 




<2ikiiii/)-i/')iiiiV'ii + 2iiuiiiii/')ii 


L't 




leading to the alternative bound flT^ . 

To show (1201) . let t G [0,T] be hxed. Suppose ■0 7 ^ 0. Then for any e > 0, we may 
choose /' G & such that 


Since T is hnite, for ko G sufficiently large we can choose 0 < < cxo and G 

for j = 1, 2,..., such that 




with Pj' = QpJjf Finally, from Corollary [H since L[ is hnite we can hnd ki G 
sufficiently large, with ki > max{p, ko}, such that 


((( 7 ,-[/,“>)■ « 0 |/'),^;) 


e 


2 



for all k > ki- From the above choices and taking square roots on both sides of ffTSj) . we 
have that 



< e, 


for all k > ki- Since the theorem statement holds trivially when ip = 0, this completes the 
proof. □ 


A discussion of the error bounds presented in Theorem [T] is now in order, starting with 
flT^ . The error bound on the right hand side of (1181) is the sum of three terms. The hrst 
term is a bound on the error committed by approximating e(/) by e{f') for some simple 
function f. The second term bounds the error in approximating Upip by a hnite sum of 
terms in given in ipp Finally, the last term gives an upper bound for the magnitude 
of the inner product between the error term (Ut — ® e(/') and ipp The hrst and 
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final terms are computable. However, the second term is difficult. We note that 

1/2 


u;-‘f - E 'i-’if 

i=i 


L't Li 

u:,p - 

i=i i=i 




J 

Li 




i=i 


2x 1/2 




(fo9j,t,o)rp(fl9j^t,l) 
t2-tl 


rp(h9j,t,l) I \ I I 

<Lj,t) r + 






i=i 


2\ 1/2 


( 21 ) 


where in the last line, we again use the identity flT3l) . However, this is difficult to compute 
as it involves the semigroup which acts on the inhnite dimensional space T-Lq. Thus 

to alleviate this difficulty we now turn to the alternative bound ([T9]). 

The hrst and third terms of (IT^ are the same as for (fTSl) . However, for the second 
term we have the identity 




i=i 







2iil|e(g',,)||S 


••• 


ti+i—tt 




1/2 




( 22 ) 


derived in the same manner as (ED. However, unlike fl2T|) . fl22|) involves only the semigroup 
Q ]2 the hnite dimensional Hilbert space Thus it is a quantity 

that will be much easier to compute. All that remains now is to construct a suitable 
approximation //''■ ^ to 'ifj. One way to do this is to choose i/i' ^ to locally minimize the 
right hand side of fl2^ or, equivalently, the term under the square root. Unfortunately, 
although this is in principle possible, it is in general a challenging and computationally 
intensive non-convex optimization problem. We will demonstrate this optimization in an 
example that will follow. 


3.3 Subspace truncation examples 

3.3.1 Kerr-nonlinear optical cavity: 

Consider a single-mode Kerr-nonlinear optical cavity coupled to a single external coherent 
held (m = 1), which is used in the construction of the photonic logic gates presented in |12] . 
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Let 'Ho = (the space of infinite complex-valued sequences with \^n\‘^ < oo) which 

has an orthonormal Fock state basis {|n)}„>o. On this basis {|n)}„>o, the annihilation, 
creation, and number operators of the cavity oscillator can be dehned (see, e.g., [T6] ) 
satisfying 

a \n) = y/n |n — 1), a* \n) = -s/n +1 |n -|- 1), a* a \n) = n\n), 

respectively. Similar to Examples 14-15 of |T6], we set "Do = span{|n) | n G The 

Kerr-nonlinear optical cavity can be described by 

S = I, L = '/Xa, H = Aa*a + xa*a*aa. 


where A, A, x > 0. We now show that Conditions!] [a) and[ ^a) hold for the Kerr-nonlinear 
cavity. 

Consider H*^^^ = span {|n) | n = 0,1, 2,..., A;} and a system approximation of the form 




= Pa,(k)HP 




(23) 


Conditions!] 1(6) and !^ 6) [hold immediately because is hnite dimensional. Note 




that H^^^ C Vq. Hence, Assumption!T]holds. Recall that = ran ( 

From the {S,L,H), we see that = span{|fc-|- 1)}. Now consider = span{|A;)} 

and = I (A/c + \a — P\‘^). We will now show that Assumptions !2| and !3] hold for the 
Kerr-nonlinear optical cavity and the approximation. 


)■ 


Assumption !^fa)! Note that 

(£'“■")-£(''^“'’))Ik<«=/3VAo‘Ik(.). 

holds with = v'A(/fc+ 1)|^|. 


Thus, we have that Assumption C (a) 


Assumption !m(6)[ This assumption follows for the dehned because 
ker ( |^(fc)) = span {|n) | n = 0,1, 2,..., fc — 1} for all a, /? G C. 


Assumption !Mc)[ First note that PjQ(k) = 0, P^{k)L^^'^P^(k) = 0, and 
P^(k) = 0. Also, for any u G H^^), Pf^wP^(k)U = 

XkP^{k)U. From these identities and the fact that is self-adjoint, we have for any 
u that 

= {Xk + |a|2 + |/3|2 - 23fJ{a*/?}) \\P^ik)uf 

-l- 3? (^'\/X(3PjQ(k)a P'j^(k)Qii^{k)U, Pit^{k)U^ 

= —g{k, a, /9) ||P^(fc)M||^ -|- h{k, a, (3, u). 
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Notice that \a\^ + |/3p — 23fJ{a*/9} = \a — /3p. Also, note that 


^ 1 P-uWqicWU, Pj^(k)U 


< l/?l^/A|| Pjt^(k)a* Py^(k)QK^(k)U\\ ||P^(fc)M 

< \(3\V^\\u\\ \\P!cik)U\\ . 


Therefore, we have that Assumption [ ^c)| holds for the dehned /C^^^^and the defined 7 ^“^^ 
with = IPly/Xk. 


Assumption For any u E Pq, HPj^(k)U = (A(/c + 1) + xk{k + 1)) Pj^(k)U. This 

implies that Pj^(k)u)} = 0. Also, L*LPj^(k)U = \{k + l)Pji^(k)U 

for any u E 'Hq- Therefore, we have that 


rnio-P) p 


= “2 “ 2SRK/3} 

+ 3? I {^a* - a*a) PmWU, 

2 

= —g{k,a,l3) Pj^(k)U + h(t, A;, a,/?, n). 

Similarly to the previous derivation, using that is a contraction, we see that Assump¬ 
tion holds for the dehned /C^^^and the dehned 7 ^“^^ with = \a\^/X(kX^ + 

WWKk + 2). 


Assumption [3t From the dehned 'y^^\ we see that this assump¬ 

tion holds for any r, s G Z+ such that r -f- s > 3. 

Finally, because Assumptions [T}13] hold. Lemma IH Corollary [T], and Theorem [T] can be 
applied to obtain error bounds on the hnite-dimensional approximations. 


Numerical example of the Kerr-nonlinear optical cavity: Consider A = 25, A = 

50, and y = —A/60 (these parameters are used in [12] )• We also consider an input held 
with a constant amplitude of a = 0.1 for t E [0,T] with T = 5. Let -0 = |0) 0 |al[o,t]). We 


will now compute an error bound on 


Ut-U, 


(k) 




for diherent values of k. From flTS]) 


(established in Theorem [T|) and the fact that al[o,i] G ©' is a simple function, we have that 


Ut-U, 


(k) 




< 4 






2 E E 

i=i *=i 
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Table 1: Numerical computation of error bounds on 
nonlinear optical cavity 


n-t/f’l i' 


for the Kerr- 


k 

Error Bound 

19 

0.2366 

29 

0.2115 

39 

0.1970 

49 

0.1872 

59 

0.1799 

69 

0.1742 

79 

0.1696 

89 

0.1658 

99 

0.1625 


where /' = al[o,t], = Ej=i V'yt, and ^ = u'j ^ ® e{g'j f) with g'^ ^ G &. Now, using 
we have 

2 


uri’--4’', =IIV>IP + IWf 

L't 

- 25 ;i|e(g',)||s {(0 

t = l 


(fc;a,^'(0)) (fc;a,^'(l)) 


rp\ ■> 




t2—tl 


rp'. ir-j 








(24) 


where /9j(*) = for any s G To find an appropriate 'ip^ we set up the 

cost function 

1/2 


Lt 


AW)= |l + IWf-2^l|e(g/,)l|S|(0 

i=i 


(A:;a,^'(0)) (A:;a,/3'(1)) 


to 




t2-tl 






t-ti 


u 




(25) 


for any ® then choose and for j = 1, 2,..., L( such that 

it is a local minimizer of Jk- 

For computational simplicity, let us fix = 1 and take tj+i — ti = 0.5 for all i = 
0,1,2, With t = ti^i = 5, we then have that i = 9. We take an initialize 

guess at = |0) (8) e(al[o,T])- Using the general purpose unconstrained optimiza¬ 
tion function fminunc in Matlab, a local minimizer 'ip[ = was found as 

u'pt = 0.9999 |0) - (0.0024 - 0.0094?) |1) - 0.000112) and g[ ^ = (0.0866 + 0.0462?)l[to,n) + 
(0.0882 -I- 0.0471?)l[t^_T)- The corresponding cost is T('0)) = 0.0096. Consider r = s = 2. 


Using error bounds on 
Recall that the dimension o 


the reduced su 


for various values of k are shown in Table [H 
□space is A; -|- 1. 


3.3.2 Atom-cavity model: 

Consider a three-level atom coupled to an optical cavity, which itself is coupled to a single 
external coherent field (m = 1) [TI]. Let "Ho = <8) We will use |e) = (1,0,0)''', 
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1+) = (0,1, 0)"'", and |—) = (0, 0,1)^ to denote the canonical basis vectors in C™. We also 
consider |n), denoting the normalized n-photon Fock state of the system, as the basis vector 
of (as in the previous example). Let Pq = ® span{|n) | n G Z+}. This atom-cavity 

system is then described by the following parameters: 


S = 1^1, L = I (S)VXq 


H = IX (cx+ (g) a — (T_ (g) a*), 


where A,x > 0, < 7 + = |e) (-f|, and (j_ = |- 1 -) (e| 
Condition!] 1(a) and[^ 


From Lemma 12 of [16], we have that 
holds for this atom-cavity model. 

Consider (g) span {|n) | n = 0 , 1 , 2 ,..., fc} and a system approximation of 

the form fl2^ . That is, we are approximating the dynamics of the harmonic oscillator. 

hold immediately because has hnite 

Hence, Assumption [ 1 ] holds. 


Recall that Conditions []|( 6 ) and 
dimension. 


We also see that C Vq. Hence, Assumption [T] holds. Similar to the 
previous example, we note that <g) |/c -|- 1). Now consider <g) |/c) and 

7 ^“^^ = A (A/c + \a — /3p). Using similar derivation to the previous example, we will now 
show that Assumptions [ 2 ] and [3] hold for the atom-cavity model and the approximation. 


Assumption [^a)[ Note that 

(£'“'=) - \ar (/3'^ + -w-) ■ 

Thus, we have that Assumption holds with = ^/k + 1 ^|/5|\/A -|- y 


Assumption [^&)[ This assumption follows for the dehned because 

ker (|^(fc)) = <8 span {|n) | n = 0,1, 2,..., fc — 1} for all a, /9 G C. 


Assumption [ 3](c)[ Note that = 0, = 0, 

0 /c('=) — 0, and P^(fc)iLP^(fc)Q;(;(fc) — *xP»c(fc) 0 '- 0 *P^(fc)e/c(fc)- Also, for any 

u G PfQ(k)L^^"> L^^^PiQ(k)U = \kPfQ(k)U. From these identities and the fact that is 
self-adjoint, we have for any u that 

(P^(.) + P'H(^)eK.w)u, P^(k)u)] 

= (Ak + |a |2 + |/?|2 - gfj{W/3}) IIP^wmII' 

-|- 9? I (^P^^ik) ^/3\/A -|- a*P-^(k)QK^{k)U, PfQ{k)ii^ I 

= —g{k, a, (3) ||Py^(fe)M||^ - 1 - h{k, a, (3, u). 


Noticing that lap -|- |/3p — 23fJ{a*/9} = (a, a) - 1 - (/9, {3) — (a, {3) — (a, (3) 
note that 


3? (P 


< 


^/c('=) Cb*P'^(k)QK2(k)U, PK^(k)U^ 

PjcW (/^v^ + xa_ ) a*P^(k)Q^(k)U ||P^(fc)'a|| 


<Vk(x+ l|w|| llPcw^ll • 


a — /3p. Also, 
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Therefore, we have that Assumption [ ^c)| holds for the defined and the dehned 
with = '/k(^x+ l/^IV^)- 


Assumption For any u G "Ho, note that L*LPj^{k)U = X{k + l)Pj^(k)U and 

iHPj^(k)U = i (A(/c + 1) + xk{k + 1)) Pj^(k)U. Also, Thus, we have that 

= + 1) + |a|2 + |/3|2 - 23fJ{a*/3}^ Pj^ik)U 


+ 3? I {/3a* - a*a) Pmwu, 

+ 3? I (cr_a* - a+a) Pj^.k^u, Pj^,k)U 


= -9{k,(y,P) 




+ h{t, k, a, f3, u) 


Similar to the previous derivation, using that is a contraction, we see that As¬ 
sumption holds for the dehned and with = yJk+1 (^x+ |Q;|V^j + 

Vk + 2 (^x + 1 / 51 • 


Assumption [S From the dehned we see that this assump¬ 

tion holds for any r, s G such that r -t- s > 3. 

Finally, because Assumptions [Tll3] hold. Lemma IH Corollary [1], and Theorem [1] can be 
applied to obtain error bounds on the hnite-dimensional approximations. 


4 Error bounds for adiabatic elimination approxima¬ 
tions 


When an open quantum system comprises subsystems evolving at two well-separated 
timescales, the system dynamics can be approximated by eliminating the fast variables 
from the model description. This method is known as adiabatic elimination in the physics 
literature and singular perturbation in the applied mathematics literature. In this section, 
we will establish error bounds for this type of hnite dimensional approximation of open 
quantum systems, when the slow subsystem lives on a hnite-dimensional subspace. 


Let satisfying ([2]) describe the time evolution of the original Markov open quantum 
system to be approximated. We set = "H and = V. We also let Ut satisfying 
(fTj) describe the time evolution of the adiabatic elimination approximation dehned on a 
hnite-dimensional subspace T-Lq C "H. In this setting, Vq = Hq. We note that Conditions 

I® 


and 


hold immediately because 'Hq is hnite dimensional. As in 


we assume 


the following. 
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Assumption 4 (Singular scaling). For j = 1, 2,..., m, there exists operators Y, Y*, A, 
A*, B, B*, Fj, F*, Gj, G*, Wij, W*^ with the common invariant domain V such that 

- if' 

2 / ^ i ^ 
i=l 

Assumption 5 (Structural requirements). The subspace TLq C FL is a closed subspace such 
that 

1. Vo = Pno'D C V 

2. Y = 0 on P 

3. There exists Y, Y* with the common invariant domain V so that YY = YY = P^±. 
4- FjP-Ho = 0 on V for all j = 1,2,... ,m 

5. P-h^AP-Uo = 0 on V. 

Assumption 6 (Limit coefficients). The approximating system operators {S,L,F[) are 
such that 


= eY + kA + B, Lf'>* = kF,+Gj, Sf* = Wij. 


m 

+ *j) 

£=1 

L- = {g, - AYF,) 

H = ^ (b - AY A) Pnl] . 


Let us now re-state an important result. 

Lemma 5 ( [T6l Lemma 10]). Suppose Assumptions\^l^ and\Bhold. The linear operators 
S, L, F[, defined in Assumptions^ have the common invariant domain Vq. The operator S 
is unitary and H is self-adjoint. 


For any a, (3 E C™ and any n G P, let us define 

( m m 

/ I |2 I I /0|2 ^ 

+B + Y,G,!},+ Y^ atW,, (ft - G*) 




V. 


i=i 




We now introduce an additional assumption required in obtaining our error bound results. 
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Assumption 7 (Boundedness of operators). For any k G Z_|_ and any a,/? G C™', we have 
that 


Ml 


{k\OL^^) _ 


M. 


(fc;Q,/3) 


Pu^ - - (^5^/3) _ j 

i j £(«/3) 

+ _|_ ^y4(“/3) _ YP-h^ 


Wo 


are /imfe (^re., 0 < < cxd^. 


Lemma 6. Suppose Assumptions\^ 0 and\B hold. Then for any k G Z+, an?/ a,(3 ^ C™, 
any u E Hq, and any t > 0, it holds that 


rjn(^k]<y.j3'j rjn 


(a/3) 


U 


< ^ ( 2Mf+ tM^ 

fC 


(k-,a,IS) 


\U\ 


(26) 


Moreover, andT^^’°‘^^ are also norm continuous for each a, (3, with || (/—T/“^^)n|| < 

A|“^^(t)||n|| and ||(/— < N 2 ^'°‘^\t)\\u\\ for some continuous nonnegative func¬ 
tions and then 


F 


{k;al3) 


(a/S) 


-F 


(27) 


-.t ju <^(Mf’“’^)(Aj“^^(t) + Af“^)(t))+tMf“’^^)||^ 

Proof. First note, from the definition of a strongly continuous semigroup, that 

1. T(S“^^ = = I for all fc G Z+, 

2. ^ C^°‘^)u for all u eTLq (since TLq is hnite-dimensional), 

3. |Tf for all u G Dom(£(^’"A). 

From the above properties and Assumption [5l we can write for all u G T-Lq and all f > 0 
that 


d 

dt 


F, 


(k]OL^) 


-F, 


(a/3) 


u 


_ ^£{k-,Oifi)rjn{k\Clj3) _ ^(^afi) rjn(al3) 

^ £(fc;a/3) j'y^(fc;a/3) _ y^(a/3)^ ^ ^ (£(fc;a/3) _ ^(a/3)^ j.{aP)^ ^ 28 ) 

with ( Tq“^^ — ti = 0. Since CF°^h) jg ^ generator of a semigroup on a Hilbert space. 


j^{k,ai3) _ rjn{Q-y)^^ ^ ^ Dom(£(^’"^)) (due to Assumption E]), and ( 

C^([0,t];'H) (the class of continuously differentiable functions from [0,f] to P), a unique 
solution of fl28|) exists and is given by [TH Thm 3.1.3] 

rt 


F, 


{k\a.0) 


-F, 


(a/3) 


u = I _ £(«/ 3 )) 


(29) 
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for alH > 0 and all u G "Ho- 

From Assumption 01 we now note that 




for any u E T> and any a,/3 E C™'. Using Assumptions 0] and El it has been shown 
that [ini p. 3146] 




for any u eVq and any a, [3 E C^. 

From the above identities and Assumption 01 for any v E Po, we have that [TB] 


^£(fc;o/3) _ £ia/S)'^ y 

_ I £(fc;a/3) _ £(o/ 3) _ ^£(A:;o^)y^(a/3) _|_ ^£(fc;o/3)y^(Q;/3) 

\ k k 

_ ^£ik-,a3)YP J_ ( _ y4(“/3)y^(«/3)'\ 

A:2 ^0 V J 



We also note, using integration by parts (in a similar manner to [T9| p. 898] and [20| Eq. 
(2.2)]), that 


f ji(fc;«/3)£(fc;a/3)y^(a/3)2n(o/3)^^.^ _ /^2n(fc;«^)y^(Q:/3) _ 

/ i *7” 7” \ £ ij j 

Jo ^ ^ 



where we have used that -f and -^Y= Y 


[T8] . Similarly, we have that 



^rj.(k-,afi)yp^^ l^p(a/S) _ ^(a/3)y^(a/3) j _ _ ^(a/3)y^(a/3)^ pM)^ ^ 


i 
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Now, because G 'Ho for all t > 0, by substituting the above three identities into 

fl2^ . we have that 




{k-,aj3) _ rp{al3) 
-‘-t 


U 


+ ^ “ \ 

-g y ^^(oiP)y _|_ YP-^± — y4(“/^)yy4("/5)^^ 

xTj“^)Mdr] . 

Here, we have used Assumption O that P-^^AP-^^ = 0. We then have that 




(k-,aP) rpiaP) 


1 

< - 
- k . 


Ml 


(k-,aP) 




{k\(y.f3) 


+ 


i(a/3) 


+M, 


(k-,afi) 


z 


(k-,al3) 


z 


A A 


dr 


\u\ 


The result fl2BD then follows from the fact that and are contraction semi- 


fand Tf 

groups. Moreover, when and are also norm continuous and satisfy ||(/ — 

< A^j'“^^(f)||n|| and ||(/— < A' 2 ^’“^^(f)||M|| as stipulated in the lemma, 

we have 


z: 


(fc;o/3) rpial^) 


-t 


U 


(J - Tfi 




+ / 


+ J _|_ ^y4("/3) _ ZP-^± (^pAP) — y4(“/3)yy^(«/3) 

xT^“^)Mdr] . 

From this and the contractivity of and it follows that 


T. 


(fc;a/3) _ 2 ,(o/3) 


< 


k 




(fc;Q^) / 


f-{k-,a0). 


\U\ 


This establishes the lemma statement. 


□ 
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Recall that & C L^([0, T]; C"*) is the set of all simple functions in L^([0, T]; C™), which 
is dense in T]; C™'). Similar to the previous section, let Uq = {u<^ e(/) | u E Hq, f E 

©'} and = {m 0 e(/) | u E f E ©'} for k > 0. 

Lemma 7. Suppose Assumptions\^\^\Bi an <*0 hold. Then for any k E any "01) V ’2 £ 
Uq, and any t E [0,T] with 0 < T < oo, we have that 




< 


1 ^ 


2Ml 


{k;ai{i)a2{i)) 


+ {U+1 — ti)M2 


(k;ai{i)a2{i)) 


ll^l|||l^2 


i=0 


(30) 


where to = 0, = t, and fj = Z)i=o i Moreover, if and 

j^{k,ap) norm continuous for each a, ft, with ||(/ — and 

11(1 — < iV 2 ^’“^^(f)||M|| for some continuous nonnegative functions and 

N^,oiy) ^ then 




V'1,'02 


1 ^ 
j=0 ^ 


(i)) 


(31) 


Proof. This proof follows similar arguments to the proof of Lemma H] using the result fl26|) 
established in Lemma |6l □ 


Corollary 2. Suppose Assumptions^\^\^ and^hold. For any t E [0,T] withO <T < oo, 
any ^fl = ui® e(/i), -^2 = U 2 ® e(/ 2 ) EPLo® iP, we have that 


(^(uP - uf ipui’-f < 2(l|Milllle(/i) - e(/0|t||i/>2ll + l|M2lll|e(/2) - e(/) 




i=0 


Ak-,a\{i)a2{i)) 


+ (t 


i+l 




IKIIII^2l|- 

(32) 


for any '0' = Uj 0 e(/j) E U.'P with for some i E 1^+ and a sequence to = t) < p < ... < 
ti < te+i = t such that f[ = Z)i=o /2 = Z)i=o “2(*)l[n,n+i) ■ addition, 

andTj;’^'°‘^^ are also norm continuous for each a, ft, with || (/—Tj"^^)n|| < iV^^"^^(f)||M|| 
and 11(1 —< iV 2 ^’“^^(t)||M|| for some continuous nonnegative functions and 

then 


-Ut'j < 2(^||ni||||e(/i) - e(/()||||?/i 2 || + ||m 2 ||||e(/ 2 ) - e^f^ 




pfc;cn (1)0:2 (*)) ( /\r("(*)/ 3 (h) 


(ivr 


’ (*i+i - u) 


i=0 


+ivf ;“(•)«“)) (ii+1 - i.))+(«i+i - ) ii,a; 


r(A;;ai(i)o 2 (*)) 


(33) 
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Moreover, it holds that 


lim 

k^oo 




0 . 


(34) 


Proof. This proof follows similar arguments to the proof of Corollary [T] using fl30|) and fl3T]) 
established in Lemma [71 □ 


Theorem 2. Suppose Assumptions^\^\^ and^hold. Lett) <T < oo. For any t G [0,T], 

eonsider any L[ G , any V’yo with V’j i = ^ ® ^ 7^ 0- 

Also, consider any f G &. Let £ be a positive integer and to < ti <...< ti < = t 

be a sequence such that f = Y!i=o^i'^)Mu,u+L 9'j,t = Let u G 

with ||n|| = 1, and |/) = e(/)/||e(/)|| G T (i.e., |/) is a coherent state with amplitude f), 
and = u \f). Then, 


(t/f - £/,)■ V'f < 4(111/) - i/')ii + \\u,y, - v-:i 


P e 




j=i i=i 




+ (L+i ~ ti)M.^ 


(k\a{i)yj,t(i)) 


\Witl (36) 


If in addition, and are also norm continuous for each a, 13, with ||(J — 

and || (/— for some continuous nonneg¬ 
ative functions and , then 


(uP - u,yi,f < 4(111/) - 1 /')|| + Iin'v- - /.'i 


P I 


iEEK 


j=l i=l 
r(k-,a{i)l3ii)) 


Ak-Mi)3j,t{i)) ( Ar(a(*)^h)) 








Wj.tl (36) 


Moreover, it holds that 


lim 

k^oo 



0 . 


(37) 


for any t G [0, T] with 0 < T < cx). 


Remark 2. As with Theorem [21 a stronger result of strong convergence uniformly over 


compact time intervals hmfc_).oo supo<t<T 


- Ut 

Theorem 11] for adiabatic elimination based on a Trotter-Kato theorem, but without error 
bounds for finite values of k. 


fj 


= 0 has been established in 


Proof. This proof follows similar arguments to the proof of Theorem [T] using f[32|) and f[3^ 
established in Corollary [51 □ 
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4.1 Adiabatic elimination examples 


4.1.1 Elimination of a harmonic oscillator: 

Consider a class of open qnantnm systems that comprises a finite-dimensional atomic 
system conpled to a harmonic oscillator which is driven by m external coherent fields 
(originally presented in [lEI)- Let H ^ where 'H' is a finite-dimensional Hilbert 
space and is the space of infinite complex-valned seqnences with Similar 

to Example 13.3.11 let {|n)}„>o be an orthonormal Fock state basis of On this basis, the 
annihilation, creation, and nnmber operators can be defined (see, e.g., [16]) satisfying 

a \n) = \/n |n — 1), a* \n) = \/n +1 |n -|- 1), a* a \n) = n\n), 

respectively. Following [16], we choose the dense domain T) = 'H'®span{|?7,) | n G Z+j. Let 
ns define | Consider the system operators 

defined snch that 


j^ik) ^ ^ ^ 0 a* -1- Eoi 0 a) -t- Eoo 0 I, 

= kEj 0 a* -f Gj 0 /, 

S-f = Wij 0 I 

where Eu, Eiq, Eqi, Eqq, Ej, Gj, Wtj are bonnded operators on the finite-dimensional space 

W. 

Now consider "Hq = = "H' 0 C |0). That is, the harmonic oscillator is eliminated 

from the model as it is forced into its gronnd state (i.e., |0)) in the limit as fc —)■ oo. In 
qnantnm optics, this process is the adiabatic elimination of an optical cavity in the strong 
damping limit. Now consider an approximation system with the operators (S', L, H) which 
are defined snch that 


m 

{FliEur^E, + 6ej) 0 I. 

i=i 

L* = {G,-Eoi{Eu)-^E,)®I, 

H {Eoo - Foi(En)-'Eio} 0/. 


Here, we stress that S is nnitary and H is self-adjoint [161 Proposition 14]. 

It has been shown in [16] that Conditions [1] and [2] hold for the above systems. The 
original system satisfies Assnmption H] with Y = Eu 0 a*a, A = Eiq 0 a* -1- Eqi 0 a and 
B = Eqq 0 J. Snppose Eu has a bonnded inverse, then Assnmption E] is satisfied with 
Y which is defined snch that Y'lp ^ \n) = -{Eii)~^il> 0 In) for n > 1 and '0 £ ■ The 

operators {S,L,H) also satisfy Assnmption O It now remains to show that Assnmption [7] 
holds. 


0 . 


Note that = 

+ E,„ + Y.l.iGA + E 


Now let = Eio + E7=i Fjl^j and let = 

Tj=i^l^i 3 ~ ^*j) ■ Note that and pGP) 
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are 






bounded operators on the finite-dimensional Hilbert space l-L' for any a, /3 G C”* since i?io, 
EoO) Fj) Gj and Wij are bounded operator on Ti'. For any 'ip E Ti', we have that 

® |0) = (En)-iQ(^V ® |1), 

{YP^^AGP)yaGP))Pj ® | 0 ) = ® | 2 ), 

V ^ 

{BGP)yAGF)^ (g) | 0 ) = (g) | 1 ), 

| 0 ) = ^qG)(^E u)-^Q^)p, (^ \^) 

^ m \ 

Eoi - a:W,,F; (En)-'Q(^V ® | 1 ), 

*,j=i / 

{BGP^YP^^AGP)yaGP))iIj (g | 0 ) = ® | 2 ). 

V ^ 

From these identities, we see that Assumption [7] holds because En has a bounded in¬ 
verse, Poi, id-, fhii are bounded operators, and pAA jg defined on the finite¬ 

dimensional subspace "Ho (he., it is also a bounded operator). Thus, the conditions of 
Lemma [ 6 l Lemma [3, Corollary [2l and Theorem [2] have been verified. 

4.1.2 Atom-cavity model | |16L Example 15]: 

Consider a system consisting of a three-level atom coupled to an optical cavity. The 
cavity and the uncoupled leg of the atom is driven by an external coherent field (m = 1 ). 
Let l-i = PP" ® (A. As in the previous example, we consider the orthonormal Fock state 
basis {|n)}„>o of f'. We will use |e) = (1,0,0)"'', |- 1 -) = (0,1,0)"'", and |—) = (0,0,1)"'' 
to denote the canonical basis vectors in C^. In the basis {|e), |- 1 -), |—)}, let us define 
= |e) (-|-| and = |e) (—|. We also define and P_ = |—) (—|. Here, 

P = (gspan{|n) | n G Z+}. Under the rotating wave approximation and in the rotating 
frame of reference, the system is described by the following operators: [I 6 ] 


= J, 

= J (g) ky/^a, 


a. 


(+) 


(g) a* j -|- ik ^cr(,_ ^ 


a 


a[_ 


where 7 , 5 ^ > 0. Here, a G C is the amplitude of the external coherent field driving the 
cavity and the uncoupled leg of the atom. 

Now consider "Ho = Po = span{|-(-) 0 |0), |—) 0 |0)}. That is, the cavity oscillator and 
the excited state of the atom (i.e., |e)) is eliminated from the model in the limit as fc —)■ 00 . 
Consider an approximating system described the operators (5, L, H) which are defined as 


S={I- 2P_) 0 /, 


L = 


. 2727 +)^-) 




P = 0. 
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It can be easily verified that S is unitary and H is self-adjoint. 

Again, it has been shown in [16] the above systems satisfy Conditions [T] and [ 2 l We then 
see that Assumption 0] with Y = —'^®a*a+g ® a* — (g) aj, A = 

I, B = 0, F = ^/7 ® a*, G = 0, and W = I. Let us define 

Kj = span{|+) 0 \j), |-) 0 \j), \e) 0 |j - 1)}, for j e Z+. 


Assumption [5] holds with Y which is defined, with respect to the basis {|- 1 -) 0 \j) , |—) 0 
\j ), |e) 0 \j - 1 )}, as [IS] 


Y 


'Hj - 



lii - 1 ) 0 ■ 

0 ^ 0 , 

-9\f3 0 f . 


dj 


j{j - 1)7^ 
4 


+ j 9^- 


It can also be seen that Assumption [ 6 ] holds for the defined (5, L, H). Thus, it remains to 
show that Assumption [7] holds. Note that thus, = g. 

Then note that 


f C Pi, 

C Pa, 

p(“/^)yA(“/^)p^„ C Pi, 

C Ps © Pi, 

pM)yp^^^^M)y^(a/3)p^^ C Pa. 

From the above relations and the fact is defined on the finite-dimensional subspace 

Po (i.e., it is a bounded operator), we have that Assumption [7| holds. 


Numerical example of the atom-cavity model: Consider 7 = 25, g = 5, and t G 


[0,T] with T = 1. Let = (|—) 
the error bound on (Ut — P 
Theorem [2|) and the 


(k) 


V’ 


0)) 0 |al[o,r]) with a = 0.1. We will now compute 
for different values of k. From (1351) (established in 


act that al[o,T] G ©' is a simple function, we have that 




(fc) 






j=l 1=1 


2Ml 




+ (L+l ~ ti)Mr^ 




Wr 




where V’t = Xlpi V’p, ^md ^ 0 e{gj ^) with g'j ^ G &. Similar to Example 13.3.11 

II~ V’tll is bounded by fl2T|) . Again, to find an appropriate let us define the cost 
function J('0() as in fl25l) for any = Yl!jLi u'j^t ® ^i9'j,t) ^ We then find ip'^ that is a 
local minimizer of J by numerical optimization. 
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for the adiabatic 


Table 2: Numerical computation of error bounds on 
elimination approximation in the atom-cavity example 


k 

Error bound 

10 ^ 

0.9347 

10 ^ 

0.2957 

10 ® 

0.0309 

10 ^ 

0.0131 

10 ® 

0.0096 


{Ut - f/f) 




For computational simplicity, we £x = 5 and set tj+i — ti = 10“^ for all i = 
0,1,2, With t = T, we then have that i = 999. Note that simultaneously op¬ 

timizing over 1000 time intervals is computationally intensive. Thus, we simplify the 
computation further by optimizing sequentially over blocks of 10 time intervals at a time. 
Thus optimization is done over 100 blocks. Optimization of the hrst block is initialized 
with = (|—) ® |0)) ® e(al[o,tio]) for all j = 1,2,... ,Lt. The optimization result of 
each block is then used to initialize the optimization of the next block in the sequence. 
A local minimizer -0'^^ = ^ ® e(( 7 'J, for j = 1,2,Lt, was found using Matlab as 

before, given by = (0.1008 + 6.15492)(|+) 0 |0)) + (0.2311 + 0.1012z)(|-) 0 |0)), 
= (0.0949 - 0.d031z)(|-F) ® |0)) -F (0.1500 + 0.02702)(|-) 0 |0)), = (-0.0245 + 

0020*)(|+)(g)|0)) + (0.2012+0.031b)(|-)(g)|0)), = (0.3266-0.143k)(|+)®|0)) + (0.2928- 

0.0995*)(|-) O |0)), u'^^ = -(0.4956 + 0.010b)(|+) ® |0)) + (0.1232 - 0.059h)(|-) 0 |0)). 
The overall resulting optimized cost is = 0.0046. Using the error bounds on 


(Ut - Uf ^ 


* 


V’ 


using fl3^ for various values of k are shown in Table [2l 


5 Conclusion 

This work has developed a framework for developing error bounds for hnite dimensional ap¬ 
proximations of input-output quantum stochastic models dehned on inhnite-dimensional 
underlying Hilbert spaces, with possibly unbounded coefficients in their QSDEs. The 
framework exploits a contractive semigroup that can be associated with the QDESs. This 
gives for the hrst time error bound expressions for two types of approximations that are of¬ 
ten employed in the literature, subspace truncation and adiabatic elimination. The bounds 
are in principle computable and vanish for each t in the limit as the parameter k, repre¬ 
senting the dimension of the approximating subspace in the case of subspace truncation 
approximation or a large scaling parameter in the case of adiabatic elimination, goes to cxo. 
The theory developed was applied to some physical examples taken from the literature. 

There are several directions for further investigation along the theme initiated in this 
paper. Devising a more efficient method for computing bounds for the term ip — 

J2j'^j,t\\ for subspace truncation and \\Ut*'ip — J2j'4’j,t\\ for adiabatic elimination, beyond 
the computationally intensive optimization based approach that was considered herein, will 
be important. Tensor network methods that have recently met a lot of success for efficient 
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simulation of one dimensional many-body systems could potentially be important for this 
purpose. There also remains the question of the conservatism in the error bounds and if 
there could be tighter bounds that can be achieved by using a different set of assumptions. 
In the numerical example of adiabatic elimination, the bound fl35ll was employed rather 
than the potentially less conservative fl36ll . This is because the latter requires determining 
whether is a norm continuous semigroup and hnding the bounding function 

a non-trivial task in general that deserves further investigation. Moreover, it would be 
interesting to see if there are exactly solvable QSDE models of a physical system with an 
inhnite-dimensional system Hilbert space for initial states that he in a hnite-dimensional 
subspace of the original, against which the conservatism of error bounds can be assessed. 
The authors are currently unaware of any such exactly solvable models. 
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